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Why  study  odor  tracking? 


•  Engineer  odor  tracking 
systems 

-  Gas  leaks 

-  Hazardous  waste 

-  Explosive  devices 

-  Biological  principles 


•  Understand  biological 
principles 

-  Mate/food  finding 

-  Pest  population  control 

-  Engineered  odor  tracking 
systems 


Manduca  sexta 
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Odor  tracking  overview 


•  Odor  plume  properties 


no  strong  gradient 


-  carried  by  wind  

-invisible  •  •  •  •  . 

-turbulent  ’  •  ^  ^ 

-artificial:  hazardous  waste,  explosives 

-  natural:  pheromone,  food 

•  Tracking  agents 

-engineered:  UAVs,  UGVs,  AUVs 

-  natural:  moths,  cockroaches,  fish 
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Moth  Odor  Tracking  Studies 


Case  Western  I^seeve 

UNIVERSITY 


cameras 


EST.  1826 


pheromone 

Qniirr.P 


100  cm 


250  cm 


Traditional  2D  Anaiysis 


Odor-modulated  anemotaxis 

•  Move  upwind  (surge)  when  odor  detected,  move  downwind  (cast)  when  odor  lost 

•  Combination  of  straight  legs  and  turns  (zigzagging) 

•  Turning  controlled  by  a  timing  mechanism  (interturn  duration  of  580  ms) 

Odor  plume  altitude  is  maintained 


Bio-Inspired  Odor  Source  Localization 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


5 


3D  Analysis 


^1^ 


•  Counter-turns  vertically  with  average  ITD 
of  550  ms  (90-95%  of  horizontal  ITD) 

•  Vertical  track  width  75%  of  horizontal 
track  width 

•  Temporal  relationship  between  vertical 
and  horizontal  turns  is  unpredictable 

•  Turns  are  ambiguous  (in  both  y  and  z) 


horizontal  position  vs.  time 


o 


vertical  position  vs.  time 


Rutkowski,  Quinn,  and  Willis 
J  Comp  Phys  A,  2009 
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3D  Analysis 


Downwind  View  of  Moth  Track 

0.3r 


-0.2 


.3  -0.2  -0.1  0  0.1  0.2  0.3 

horizontal  position,  z  (meters) 


•  Turns  continuously  as 
viewed  from  downwind 

•  Turn  rate  varies 


•  Turn  direction  varies 


Rutkowski,  Quinn,  and  Willis 
J  Comp  Phys  A,  2009 
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Instead  of  using  inter-turn  timers  to  control 
horizontal  and  vertical  turning  behaviors... 

•Motion  decomposed  into  normal  (v^)  and 
tangential  (v^)  components  relative  to  wind 

•Turn  rate  of  v^,  i|j=d0/dt  depends  on  odor 
concentration 

•Magnitude  of  depends  on  odor 
concentration 


Rutkowski,  Quinn,  and  Willis 
ICRA,  2007 
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1 .  Measure  odor  concentration 

Two  linearly  responding  odor  detectors 


Rutkowski,  Quinn,  and  Willis 
ICRA,  2007 
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1 .  Measure  odor  concentration 

2.  Determine  wind  direction  (not  trivial) 

3.  Estimate  odor  plume  centerline 

Average  of  antenna  position  (in  wind- 
normal  plane)  weighted  by 
concentration  and  a  “forgetting 
factor” 
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1 .  Measure  odor  concentration 

2.  Determine  wind  direction 

3.  Estimate  odor  plume  centerline 

4.  Calculate  desired  turn  rate  itJ=d0/dt 

Turn  sharply  (6  rad/s)  if  c<c,hreshoid 
Turn  as  linear  function  of  concentration  if 

threshold  ■^^■^^saturation 

Turn  softly  (1  rad/s)  if  c>Csa,^ration 
Turn  toward  estimated  source  location 

sign(HJ)=sign(03) 
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1 .  Measure  odor  concentration 

2.  Determine  wind  direction 

3.  Estimate  odor  plume  centerline 

4.  Calculate  desired  turn  rate 

5.  Calculate  desired  normal  velocity 

Vn=30  cm/s 

0=j  i|j  dt 
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1 .  Measure  odor  concentration 

2.  Determine  wind  direction 

3.  Estimate  odor  plume  centerline 

4.  Calculate  desired  turn  rate 

5.  Calculate  desired  normal  velocity 

6.  Calculate  desired  tangential  velocity 
(Vt) 


Surge  upwind  (30  cm/s)  if  c>Cthreshoid 
Cast  slowly  downwind  (7  cm/s)  if  c<Cthreshoid 
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Real  vs.  Simulated 


Results 


What  if  odor  tracking  continues  after  reaching  source? 


Tracker  remains  in  vicinity  of  source! 
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Results 


colored  dots  -  tracker  position 

big  black  dot  -  estimated  plume  centerline 


horizontal  position,  z  (mm) 
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downwind  position,  x  (mm) 


Results 


casting  and  surging 
counter-turns  without  inter-turn  timer 


“ambiguous”  turns 
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That’s  great,  but  how  are  wind  velocity 
and  self-motion  velocity  estimated? 
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Insect  Navigation  System 


•  no  GPS 

•  no  dedicated  acceleromete 

•  eyes 

-  not  distance  estimators 

•  too  close  together 

•  fixed  focus 

•  small  region  of  overlap 

-  optic  flow 

•  head  hairs  and  antennae 

-  speed  and  direction  of  air- 
current 
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unique  values  for 

-  height  (h) 

-  groundspeed  (Vg) 

-  wind  velocity  (w) 

cannot  be  determined  from 
single  measurements  of 

-  air-current  velocity  (vj 

-  optical  flow  (A) 

However... 


height  can  be  estimated  by  assuming 
wind  is  smooth  over  short  time 
period  (~1  s) 


^  '  N 


V  '  'N  V  ‘  ^ 
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Egomotion  and  Wind  Velocity  Estimation 


Moth  flight  track  recorded  at  30  Hz 


Wind  data  from  open  field  recorded  at  30  Hz 


^  u.O|  ,  ,  ,  ,  , 

8  raw 

>  Oil,  - nnered 

V 

0  5  10  15  20  25  30  35  40  45  50 
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component  of  wind  velocity  (m 

7.  ‘  t  o  ° 

rew 

- filtered 

*  ®0  5  10  15  20  25  30  35  40  45  50 

time  (seconds) 


Simulated  airspeed  data  :  =  (Vg  -  w 

Simulated  optical  flow  data  :  A  =  (Vg  /h) 

simulated 
sensor  noise 

optical  flow:  30  Hz,  a=0.2  rad/s 
airspeed:  30  Hz,  a=100  mm/s 
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Height  Estimation  Resuits  Using 
Aero/Opticai  Fusion 


m=2,  n=50 


0  2  4  6  8  10  12 

t(s) 

Quality  of  velocity  estimation  depends  on  quality  of  height 

estimation 
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Ground  speed  Estimation  Results  Using 

Aero/Optical  Fusion 


Ground  speed  (Vg)  estimate  lags  by 
two  timesteps 

Lag  caused  by  finite  difference 
approximation  of  velocity  from 
optical  flow 
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Wind  Estimation  Results  Using 
Aero/Optical  Fusion 


Wind  estimate  also  obtained 

Lag  also  present,  but  minimal  compared 
to  noise 
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h  (m) 


Odor  Tracking  and  Aero/Optical  Fusion 


•  Fusion  of  vision  and  airspeed 
makes  odor  tracking  possible 

•  Height  estimation  fails  with  constant 
velocity  motion 


•  Motion  excitation,  as  produced  by 
the  odor  tracking  task,  ailows  for 
state  observability 


-800 


time,  t  (s) 


time,  t  (s) 
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Summary 


-  Moth  odor  tracking  behavior  in  3D  is  NOT  a  simple 
extension  of  2D  ideas 

-  Moth-like  “counter-turning”  behavior  can  be  achieved 
without  inter-turn  timers 

-  State  and  wind  estimation  can  be  performed  using  insect 
sensory  system 

-  A  “drunken  stumble”  can  actually  be  beneficial 
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BUGS  PU  \N  SUCH  CRAvZ.V 
LOOPS  AMD  ZIGZAGS.  I 


Thank  you  for  your  attention! 
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Additional  Slides 
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Odor  Plume  Model 


•  Normally  distributed  odor 
concentration 

•  Plume  narrows  with 
approach  to  odor  source 

•  Centerline  concentration 
increases  with  approach 
to  odor  source 


2.  AXIS 


Zannetti 


But  remember...  odor  plume  is  patchy 


Bio-Inspired  Odor  Source  Localization 


Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited 


31 


number  of  observations 


Odor  Plume/Sensor  Model 


Odor  concentration  at  a 
single  location 


odor  concentration 


Odor  concentration  at  a 
single  timestep 


2c(x.y.z) 
Co(x)  ' 


a  -0^  -Tz 


Randomness  simulates  turbulence 


Time-averaged  odor 


ccirtcentra 


tion 


2c(x.y,z) 


c(x.y.z 


Cn(x)  ' 


tu 


odor  s 


rce 


Normally  distributed 
concentration  profile 


wind 


odor  source 


Odor  plume  narrows  with  approach  to  the  source 
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magnitude  of  (rad/s)  sqrt(odor  concentration) 


Results 


3 


time,  t  (s) 


controls 

upwind  velocity 
turn  rate 
as  function  of 

odor  concentration 
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Spiraling  Results 


15  -30 


Turns  in  the  direction  of  the  estimated  odor  source  location 


-150  -30 


d4 
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